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ABSTRACT 

Recent Very Long Baseline Interferometer (VLBI) observations determined the distances and proper 
motions of star-forming regions in spiral arms directly. They showed that star-forming regions and 
young stars have large peculiar motions, as large as 30 km s~^ with complex structures. Such a large 
peculiar motion is incompatible with the prediction of the standard theory of quasi-stationary spiral 
arms. We use a high-resolution, self-consistent 7V-body+hydrodynamical simulation to explore how 
the spiral arms are formed and maintained, and how star-forming regions move. We found that arms 
are not quasi-stationary but transient and recurrent, as suggested in alternative theories of spiral 
structures. Because of this transient nature of the spiral arms, star-forming regions exhibit a trend of 
large and complex non-circular motions, which is qualitatively consistent with the VLBI observations. 
Owing to this large non-circular motion, a kinematically estimated gas map of our Galaxy has a large 
systematic errors of 2 — 3 kpc in the distance from the Sun. 

Subject headings: galaxies: structure — galaxies: kinematics and dynamics — galaxies: spiral — ISM: 
structure, kinematics and dynamics — Galaxy: structure — Galaxy: kinematics 
and dynamics — mascrs — method: numerical 



1. INTRODUCTION 

How the spiral arms in a disk galaxy are main- 
tained has been a long-standing mystery. Standard 
understanding is that stellar spiral arms are main- 
tained as stationa ry density waves (iLin fc Shul [1963 : 
iBertin fc LinI I1996D . In this theory, the interstellar 
medium (ISM), whose motion is slightly perturbed by 
these spiral arms, forms visible spiral arms of star- 
forming regions and young stars. Alternative theories, in 
which stellar spiral arms are transient and recurrent (e.g . 
Goldreich & Lvnden-Bell' 1965': 'Julian feToomii [l966l : 
Toomrc 1981; Sollwood & Carlborg 1984), have been also 
proposed, and there has been no conclusive observa- 
tional evidence or theoretical argument on which theory 
is right. 

Recent V ery Long Baseline I nterferometer (VLBI) ob- 
servations (|Reid et all l2009bl . and refer ences therein) 
determined the distances and proper motions of star- 
forming regions in spiral arms of our Galaxy directly. 
The distance and the three-dimensional motion of a star- 



forming region W3 OH in the Perseus spiral arm have 
been determined through the parallax measur ement us- 
ing t he Very Long Baseline Array (VLBA) (jXu et al.l 
I2006f) . Its distance fr om the Sun is 1.95±0.05 kpc and its 
peculiar motion (i.e. deviations from the circular rota- 
tion) is {U,V,W) = (17±l,-14±l,-0.8±0.5)kms-i, 
where the three components are the residual from the 
circular velocity in the direction of the galactic center 
(?7), the direction of rotation (y), and the direction per- 
pendicular to the galactic plane {W) ^ . The two compo- 
nents in the galactic plane sum up to 22 km s^^, which 
is about 10% of the circular velocity of the sol ar system, 
Oo = 220 km s"! (iKerr fc Lvnden-Belll [l98l . W3 OH 
is not an exceptional object. Table [T] and Figure [1] show 
the compilation of the results of recent high-accuracy as- 
trometric observations of distances and peculiar veloci- 
ties of star forming regions using the VLBA and VLBI 
Expl o ration of Radio Astrometry (VERA ) (IXu et al 
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Brunthaler et al.H2009t iBartkiewicz et all \2QQ^ . Most 
sources have peculiar velocities larger than 10 km s~^ 
with various directions. 

These large peculiar motions of star-forming regions 
are incompatible with the standard understanding of the 
spiral arm based on the stationary density wave theory 
(|Lin fc Shul I1964D and a standing gala ctic shock solu- 
tion in a tight - winding spiral potentia l (|Fuiimotol 119681 : 
[Robert j [l969t [Roberts fc YuanI I1970D . which predicts 
non-circular flow depending on the strength of spiral po- 

* In IXu et ahl II2006I ). the solar motion relative to the lo- 
cal standard of rest (LSR) is assumed to be {Uq,Vq,Wq) = 
(10.0,5.25,7.17) km s"! l IDehnen fc Binnevl[T99g) . and the galac- 
tic radius and circular rotation of the LSR are assumed to be 
Rq = 8.5 kpc and Oo = 220 km s^^ 
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tential, whose amplitude is ~ 8 — 12 km s ^ ev en for a 
very strong spiral (i.e. F = 5% IShu et anil973[ ). More 
importantly, the "galactic shock" solutions predict that 
the spiral perturbation can only generate a laminar flow 
(i.e. streamline flow). This picture does not match the 
observations. 

It has been recognized that ISM in external barred 
spiral galaxies, such as M51, show larg e non-circular 
motions (e.g. iGarcia-Burillo et al.l Il993f l. This has 
been also confirmed by hydrodynamic simulations of 
an is othermal gas in a bar potential (e.g. lAthanassoulal 
Il992[ ). Peculiar features of the longitude- velocity di- 
agrams of atomic and molecular ISM in our Galaxy 
are also interpreted as non -circular flows under the in- 
fluen c e of the stellar bar (iBinnev et al.' '1991|; ICombesI 
fT99lHIodri guez-Fernandez fc Combes 2008 ) . However , 
almost all previous numerical studies of gas dynamics in a 
barred-spiral potential were not realistic enough to be di- 
rectly compared with observations by the VLBI, because 
star-forming dense gas clouds were not self-consistently 
treated in their models. Either an isothermal equation 
of state with g as temperature of 10^ K was assumed 
(iFriedh fc Benj[l9M) . or the ISM was represented by 
coUis i onal particles with a ma ss spectrum tT omisaka 
[l98l IGarcia-Burillo et all 119931) . iWada fc Kodai ([200 
clearly showed how the structures of the ISM diff'er in 
the isothermal and the multi-phase gas models in a bar 
potential. 

Recently, simulations of gas disks in a flxed (i.e. time- 
independent) background stellar potential were tried 
with very high spatial resolutions, bo th with an Euler 
grid technique (jWada fc Normanll2007f ) an d a smoothed 
parti cle hydrodynamics (SPH) technique () Salt oh et al.l 
|2008() . Thanks to the high spatial and mass resolutions, 
they were able to follow the cooling of gas to quite low 
temperature (below 100 K), and thus found that gas 
disks develop complex fractal-like structures. In these 
simulations, quasi-steady distribution of gas was main- 
tained through the balance between the two competing 
processes. One is the structure formation through radia- 
tive cooling and gravitational contraction, and the other 
is the disruption through tidal interactions, heating by 
supernova (SN) explosion and shear due to a differential 
rotation. However, no clear spiral structure was formed. 
This is not surprising because the stellar disk is expressed 
by a fixed, axisymmetric potential field or by A^-body 
particles without a sufficient resolution. 

In the present paper, we combine high-resolution N- 
body and hydrodynamic calculations in order to under- 
stand the peculiar motion of star forming regions ob- 
served in our Galaxy. In Sj2l we describe our numerical 
model, in which a quasi-steady, isolated 'spiral' galaxy in 
a fixed spherical potential is produced from a featureless 
stellar and gas disk. The evolution of the stellar disk, and 
distribution and kinematics of cold gas and young stars 
are presented in SjS] This section also compares peculiar 
velocities of star- forming gas in our simulation with those 
of observations. In fjU we discuss the rotation speed of 
the star-forming gas, origins of large peculiar motions 
of star-forming gas, and how kinematic distances of the 
Galactic objects differs from the true ones. The results 
are summarised in [J5l 



2.1. Numerical Methods 

In our model, evolution of a stellar disk with the 
multi-phase gas in a halo potential, taking into account 
star formation from cold, dense gas and energy feedback 
from SNe, was solved numerically. We used our original 
A^-bo dy/gas simulation code ASURA (jSaitoh et all l2008l . 
I2009f ). Hydrodynamics was solved by the standard SPH 
methods and the smoothing length was allowed to vary 
both in space and time with the constraint that the typ- 
ical number of neighbor each particle is near A ^^b = 32. 
The artificial viscosity term (jMonaghanl Il997t ) and the 
correction te rm to avoid lar ge entropy generation in pure 
shear flows ([Balsaral |1995[ ) were used. Self-gravity of 
stars and SPH pa rticles were cal culated by the Tree with 
GRAPE method (lMa.kinolll991l) wh ich is a combination 
of the Tree method ([Barnes fc Hutl [T986[) and GRAPE 
(jSugimoto et al.lll990f ). The opening angle was set to be 
0.5. In these simulations, we adopted a software emu- 
lator of GRAPE, named Phantom-GRAPE (Nitadori et 
al. in prep.). Phantom-GRAPE is tuned by Single In- 
struction Multiple Data (SIMD) instructions which are 
equipped by modern CPUs. The leapfrog time- integrator 
was adopted. 

Metallicity-dependent radiative cooling of the gas was 
solved assuming an optically thin cooling function with 
solar metallicity that cov ered a wide range of tem pera- 
ture, 20 K through lO^K (|Wada fc Norma.nl [200l . Uni- 
form photo-electric heating of dust by the far-ultraviolet 
radiation (FUV) ob served in the solar neighborhood was 
taken into account (jWolfire et al.lll995l : iGerritsen fc Ickd 
[1991 . 

The treatment of star formation a nd the heating due t o 
the SN feedback was the same as in iSaitoh et al. r il200l . 
We adopted the single stellar population approximation , 
with the Salpeter initial mass function (|Salpetedll955l ) 
and the mass range of 0.1 — IOOM0. An SPH particle 
was replaced with a star particle following the Schmidt 
law (jSchmidtlll959[ ) with a local star formation efficiency 
(C* = 0.033) in a probabilistic manner, if criteria (1) 
riH > 100 cm-3, (2) Tg < 100 K, and (3) V • 1/ < are 
satisfied. Here we assumed the dense and cold gas as the 
potential site of star formation. Adv antages of this set o f 
star formation criteria are as follows (|Saitoh et al.ll200l : 
(1) simulations reproduce realistic structures of the ISM 
and young stars, while the set of criteria (tijj > 0.1 cm"'^ 
, Tg < 15000 K, and V ■ v < 0), adop ted in many pre- 
vious studies of galaxy formation fe.g.. iGovernato et al.l 
l2009f ) do not reproduce those structures, and (2) global 
(galactic) star formation rate is not directly proportional 
to the local star formation efficiency, C* , but largely con- 
trolled by the global evolution of the ISM from reservoir 
to dense gas where st ars are formed (for details, see §5.2 
in lSaitoh et al.ll2008f l. The second point is specially im- 
portant for numerical simulations of galaxy formation 
and evolution since our result does not depend strongly 
on the value of C* used. We implemented type-II SN 
feedback, where the energy from SNe is injected to the 
gas around the star particles in the thermal energy. Each 
SN releases 10^^ ergs of thermal energy for the surround- 
ing ISM. 



2. METHODS AND MODEL SETUP 



2.2. Initial Model 
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TABLE 1: Peculiar motions (deviation from a galactic circular rotation) of star forming regions in our Galaxy 
measured by VLBI astrometric observations of H2O and CH3OH masers obtained by the VLBA and VERA. D [kpc] 
is distance from the Sun, C7, V , and W [km s~^] are three components of peculiar motions. U is toward the Galactic 
center, V is toward Galactic rotation, and W is toward the North Galactic Pole. Galactic radius and circular rotation 
of the LSR are assumed to be Rq = 8.5 kpc and Oo = 220 km s~^. The flat rotation cu rve is assumed. The solar 
motion relative to the LSR is (C/q, VpuWp, ) = (7.5, 13.5, 6.8) km s~ ^ (jFrancis fc Andersonii2009i . hereafter FA09) and 
{Uq,Vq,Wq) = (10.00,5.25,7.17) km s^^ (|Dehnen fc Binnevlll998l . hereafter DB98). 
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Fig. 1. — : Spacial distribution and peculiar velocities of the star-forming regions in our Galaxy (sec Table 1). a) 
Peculiar velocities of star- forming regions in our Galaxy, based on the solar motion relative to the LSR of FA09. b) 
Same as panel a), but the solar motion relative to the LSR of DB98. 



The procedure to prepare the initial model is as fol- 
lows. We first place a pure TV-body stellar disk of an 
exponential profile in the center of a static dark-matter 
(DM) halo potential. Details of the halo potential and 
the disk model will be given in i i2.2.1l Then, we let the 
stellar disk evolve for 2 Gyr. After a bar and stellar spi- 
rals have developed well, we added a gas component to 
the stellar disk (§2.2.2). We used this combination of 
fixed DM potential and stellar and gas disks as the ini- 
tial condition of the simulation for the gas disk-stellar 
disk-static halo model. 



2.2.1. DM Halo + Stellar Disk Model 

We adopted a fixed, spherical DM halo as a host of an 
iV-body stellar disk. The density profile of the DM halo 
follows the Navarro-Frenk- White (NFW)-profile: 

r/rs(l + r/TsY 

where 

^ _ Ml C'nfw'^ 

47rJ?h^ + C'nfw) + Cnfw/ (1 + Cnfw) ' 

Cnfw = Rh/fs, (3) 
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C'nfWi Mil, and _Rh are the concentration parameter 
of the DM halo, the mass, and the vhial radius, re- 
spectively. The concentration parameter is set to be 
Cnfw = 5. The mass is Afh — 6.3 x 10^^ Mq and the 
virial radius is i?h = 122 kpc. The values of A^h and 
i?h arc determin ed by using the spherical collapse model 
(Equation (2) in lMo et al.lll998[ ) with a circular velocity 
of the halo Vc = 150 km and a formation redshift of 
the halo zi = 1. 

The initial condition fo r the stellar disk is generated 
using Hernquist's method (|Hernauistj|1993[ ). The density 
profile of the stellar disk is given by 



TABLE 2: Model parameters of DM halo + stellar disk 



Ms, 



exp(-i?/i?sd)scch^(z/zsd), (4) 



where Mgd is the mass of the stellar disk, R^d is the 
radial scale-length, and Zgd is the vertical scale-length, 
respectively. The vertical velocity dispersion CTz is related 
to the surface stellar density Sgd and Zgd: 



=^Gz,dSsd(i?), 



(5) 



where G is the gravitational constant. The radial velocity 
dispersion ctr is also assumed to be directly related to the 
surface density: 



a^(i?) = Aexp[~i?/i?,d] 



(6) 



The normalization constant A is fixed in such a way that 
at some reference radius i?rof(= 1.5-Rsd) the radial ran- 
dom velocities are Qrof times the critical value needed to 
stabilize a differentially rotating disk against local per- 
turbations: 



■cf- 



3.36GS«iGRref) 



(7) 



where K{R^cf) is the epicyclic frequency at i?rcf- The 
azimuthal component of the rotational velocity (V^) and 
its dispersion (cr^) are foimd using the asymmetric drift 
and the epicyclic approximations: 



Kir - <^R 



2R 

Rsd 



402 



- 1 



4 



4172 



(8) 
(9) 



where Vdr and = Vdr/R are the circular velocity and 
the angular frequency at a given R, respectively. The 
model parameters of the DM halo and the A^-body stellar 
disk are summarized in Tabled 

As suggested by many previous studi es (e.g. 
lOstriker fc Peebleslll973l : lEfstathiou eranil982f) . a mas- 
sive disk exper i ences a bar instability.According to 
lEfstathiou et al.l ()1982t ). who used A-body simulations 
to investigate bar instabilities of exponential disks cm- 
bedded in fixed DM haloes, the bar instability takes place 
in the stellar disks when the criterion 



y^GMsd/ Rsd 



< 1.1 



(10) 



is satisfied. Here Knax is the maximum rotation velocity 
of the disk. In our model, the value of eb is ~ 0.8 for 
the initial state. Thus our model is unstable against bar 
formation. 



Total mass of halo (Afi^) 
Virial radius of halo (Rh) 
Concentration parameter (Cnpw) 
Total mass of stellar disk (M^^) 
Scale length of stellar disk {Rsd) 
Scale height of stellar disk (zgd) 
Toomre's Q-value at i? = 1.5 R^d (Qrcf) 



6.3 X 10^^ Mq 
122 kpc 
5.0 

3.2 X 10^° M0 
3.0 kpc 
0.30 kpc 
1.0 



2.2.2. DM Halo + Stellar Disk + Gas Disk model 

We added a gas component to the barred spiral stellar 
disk in which global quantities such as the radial surface 
density profile, Q parameter, and velocity dispersions 
did not significantly change. The gas disk has a radi- 
ally exponential density profile, whose total mass is set 
to Mgd = 0.1 Mgd, and radial and vertical scale- lengths 
of the gas disk are set to i?gd = 2 R^d and Zgd = 0.2 kpc, 
respectively. The gas disk initially has the same circular 
velocity determined by the mass distribution. Although 
the velocity dispersion of 10 km s~^ is added for the verti- 
cal direction, there is no effect on our results. The initial 
temperature of gas is set to be 10^ K. Model parameters 
of the gas disk are summarized in Table [3l 

TABLE 3: Model parameters of the gas disk 

Total mass (Mgd) 3.2 X 10^ Mq 

Scale length (JJgd) 6.0 kpc 

Scale height {z^d) 0-2 kpc 

Initial vertical velocity dispersion 10 km 

Initial temperature (Tin) 10^ K 



2.3. Resolutions 

Total numbers of old stars and gas particles are 3 x 10^ 
(As) and 10^ (Agpn), and particle masses are ~ 1.1 x 
IO^Mq (tos) and 3.2 x 10^ Mq (rngpn), respectively. 
With this gas particle mass, gravitational fragmenta- 
tion of den se gas clumps with > XQ^Mq is resolved 
(|Saitoh et a l. 2008). 

For calculations of gravitational force, we chose a soft- 
ening length e = 10 pc. This value is small enough to 
resolve the three- dimensional structure of a disk galaxy. 
iHernguis^ (|1987[ ) showed that the softening length e 
should be smaller than the mean inter-particle distance 
i at the half mass radius, in order to achieve a good ap- 
proximation of the gravitational force calculation. For 
an exponential disk, I « 2.Q2>{N^R~^ z~^)~^/^ , we ob- 
tain £ « 35 pc as an upper limit on the softening length. 
However, although the smoothing length is allowed to 
vary time and space, the typical smoothing length is a 
few times 10 pc in our simulations. Therefore, a spacial 
resolution in our simulation is estimated as a few 10 pc. 

3. RESULTS 
TABLE 4: Resolutions 



Initial number of star particles (A^s ) 

Mass of individual star particles (mg) 

Initial number of SPH particles (A^sph) 

Initial mass of individual SPH particles (mspn) 

Gravitational softening length (e) 



3 X 10" 

1.1 X 10*^ M© 
10'' 

3.2 X 10^ Mq 
10 pc 
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3.1. Evolution of Stellar Bar and Spirals 

Figure [5] shows the evolution of the stellar disk from 
t = 1.1 - 1.5 Gyr {t = means the time when the gas 
component is added to the pre-evolved stellar disk). The 
bar structure, which is seen as two distinct enhancements 
at </) ~ 90° and 270°, R < 5 kpc, stably exists, but spi- 
ral arms in the outer disk show very complex behavior. 
The structures appeared in a co-rotating coordinate of 
the bar in the plots, we can see that most of the arms 
drift downward. They drift down faster in the outer re- 
gion than in the inner region, indicating that they move 
with the local circular speed. The pattern and number of 
stellar spiral arms change within a few rotational periods 
of the galaxy (several 100 Myr). Arm-arm merging and 
arm break-ups frequently occur. 

In order to clarify the evolution of spiral arms in detail, 
we mark stars in two trailing spiral arms at t 1.1 Gyr. 
Their trajectories are shown in Figure[2l The marked red 
stars (Gl stars) and the marked blue ones (G2 stars) ro- 
tate clockwise on the figure due to epicyclic motions, but 
their radial oscillation amplitudes are as large as ~ 2 — 3 
kpc. From t = 1.1 to 1.2 Gyr, the trailing arm marked 
with the G2 stars spreads as the stars are dispersed. The 
other arm marked with Gl stars breaks up into two arms: 
one merges with the outer arm marked with the G2 stars. 
Then, a new leading arm « 200°) forms hy t = 1.19 
Gyr. After t = 1.2 Gyr, the new leading spiral is wound- 
up {t = 1.2 — 1.3 Gyr) and grows into a single trailing arm 
at i = 1.4 Gyr. This new trailing arm consists of the Gl 
and G2 star s. This is expected in the swing am plification 
mechanism (lToomrdll981l : lAthanassoulalll984l ). 

As described above, most spiral arms in the model are 
not steady, but wound-up, and eventually dimmed. A 
spiral arm often merges with other arms, or breaks up 
into two or more spirals with coherent motion of the in- 
ternal stars. Stars in a part of the arms orbit together 
with the arm in about 1/8-1/4 rotational period at the 
radius. In other words, the behavior of stellar spirals are 
closer to 'material' arms for t he period rather than to 
quasi-stationary density waves (iLin fc Shulll96l . where 
orbital motion of each star is independent of the kine- 
matics of density waves. Interaction between arms and 
the splitting of arms are sources of 'swing amplification', 
followed by growth of new trailing spiral arms. As a 
result, spirals are no longer long-lived, but rather, are 
recurrent phenomena. Although the theoretical idea of 
swing amplification is based on a lin e ar approximation 
(e.g. iGoldreich fc Lvnden-BeTi' '1965'; 'J ulian fc Toomrd 
1 196a iToomrd Il9~ l: Athanassoula 1984h, the mecha- 
nism found here is highly non-linear, in a sense that 
the radial oscillation amplitude of the stars is not small 
compared to the system size. The recurrent spirals 
were rep orted in previous two-dimen sional TV-body sim- 
ulations ( |Sellwood fc Carlbergj [198^ . However, three- 
dimensional, high-resolution calculations are necessary 
to study long-term behavior of the recurrent spirals with- 
out being affected by numerical two-body relaxation (J. 
Baba et al in prep.; M. Fujii et al. in prep.). 

3.2. Kinematics of Young Stars and Star Forming Gas 

Figure [3] shows a snapshot of the system at t = 0.75 
Gyr. A face-on view of the stars is presented in the 
top left panel. One can see a few grand-design spirals. 



a bar with the length of about 5 kpc, and a central 
bulge-like structure. This bulge resembles the bulge of 
our G alaxy observed by the COBE /DIRBE (|Dwek et al.l 
Il995f ) in the edge-on view (top right panel). Clusters of 
young stars (bottom left panel) roughly trace the back- 
ground spirals of old stars, and the multi-arm spiral fea- 
tures of the cold gas (bottom right panel) appear. Com- 
pared with the background stellar spirals, gas arms are 
much more filamentary, but they are not hydrodynami- 
cal shocks as pointed out by recent hydrodynami c sim- 
ulatio ns in a fixed spiral potential fWada 2008 : .Dobbsl 
[2008h . Some gaseous spirals are clearly seen even in 
places where the background stellar spirals are rather 
weak. These multi-arm spiral structures of the ISM are 
morphologically similar to dust emissions in nearby spi- 
ral ga laxies recently reve aled by the Spitzer Space Tele- 
scope (|Bendo et al.ll2008D . In previous simulations of disk 
galaxy formation in which the radiative cooling below 
~ 10** K is not solved (see e.g.. iGovernato et al.ll2009t ). 
the resultant structure of the gas disk is much smoother 
than that in real galaxies. 

The left panel of Figure |4] shows peculiar velocities of 
young stars overlaid on old stars at t = 1.25 Gyr, and 
in the right panel, peculiar velocities of the high-dense 
cold gas (riH > 100 cm~^ and Tg < 100 K, hereafter star- 
forming gas) are shown on top of the gas density. Here, 
the peculiar velocity (vpec) is defined as deviation from 
the true circular rotation determined by the azimuthally- 
averaged gravitational field. Many of the star-forming 
gas clumps and young stars have peculiar velocities of 
20 — 30 km s""'^, and the directions of the arrows are ran- 
dom. An average rotation speed of the star-forming gas 
and the true rotation curve are shown in Figure O The 
average rotation speed agrees well with the tr ue rotation 
curve. Contrary to the recent statement by iReid et all 
fi2009bi . in our model, there is no clear tendency that the 
star-forming gas moves slower than the galactic rotation 
(see SH). 

3.3. Pseudo-observation of Peculiar Velocities of Star 
Forming Gas 

The left panel of Figure [6] shows an example of pseudo- 
observation of the peculiar velocities of star-forming gas 
relative to the "SUN". The position angle of the "SUN" 
is assumed to be 25° from a major axis of the stellar 
bar, and its galactocentric radius {Rq) is assumed to be 
8.0 kpc^. A flat rotation curve with the circular speed 
at the "SUN" is assumed. The complex distribution of 
large non-circular motions as seen in the observations 
(Figure [T]) are roughly reproduced. In the simulation, 
most of the star-forming gas clouds inside the solar circle 
rotate more slowly than the galactic rotation by —10 
to —20 km s~^ (represented by upward arrows), whereas 
they rotate faster (i.e. V ~ +5 to +10 km s^*) in the 
outer region. These systematic flows are caused by the 
difference in the true rotation curve from the flat one. 
Contrary to the simulation, the observational data both 
inside and outside the solar circle show slow rotations, 
especially in Figure [Hb). 

^ There are many observational suggestions that the solar system 
is located around 20 — 30° (relative to the Galaetic rotation) from 
the major axis of the bar. The value of Rq is restricted to a range 
of 7 < -Ro < 9 kpc, where the observed position-velocity diagram 
is roughly reproduced. 
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Fig. 2. — : Evolution of normalized stellar density (density contrast) in the plane of aziniuthal angle and radial distance, 
with the angular speed riframo = 25.6 km s^^ kpc~^ which is equal to the pattern speed of the stellar bar. Only the 
over-dense regions are colored and white regions have the density the below the average value. The average density 
is obtained by averaging the stellar density for an azimuthal direction at each radius. Red marks and lines represent 
trajectories of the stars (Gl stars) constituted one spiral arm at t = 1.10 Gyr, blue ones represent the stars (G2 stars) 
constituted another spiral arm also at t = 1.10 Gyr. 



More quantitative comparison between the observa- 
tions and the simulation is shown in the right panel of 
Figure [51 Black triangles indicate the star-forming gas 
shown in left panels. Open circles and crosses indicate 
the observed star-forming regions, whose peculiar veloci- 
ties are calculated by assuming the solar motions of FA09 
and DB98, respectively (see Table 1). The simulation 
data points are distributed in y/"(4,sun ^ —0.2 to -1-0.1 
and slightly biased toward V <Q. Observed data points, 
on the other hand, appear in V/Vc^sun ~ —0.15 to -1-0.05, 
and apparently they are slower than the galactic rotation 
(i.e. V <Qi). This slow rotation appears more clearly in 
crosses (DB98) than in open circles (FA09), because the 
solar motion relative to the LSR in crosses (DB98) is 
Vq = 5.25 km s"\ rather than 13.5 km s^Mn FA09. 

Figure [T^a) shows another example of the pseudo- 
observation from the "SUN" at i?o = 7 kpc, it shows no 
net rotation {V k, 0). The lengths of the arrows inside 
the solar circle become slightly short, while those outside 
the solar circle slightly elongated. In contrast, if we place 
the "SUN" further out, i?o = 9 kpc, the simulation result 
shows significant bias to F < 0. Therefore, an average 
value of V depends on i?o- In addition, comparing Fig- 
ure [6] with Figure [Tl^b), details of the distribution change 
through time through the evolution of spiral arms, but 
the width of the distribution is qualitatively unchanged. 



While the observation shows a clear tendency of inward 
radial motions (fj > 0), our simulation results in all cases 
(Figures m Elja) and[71Jb)) do not. One should note that 
the observed area does not cover the whole galactic disk. 
In particular, the data points with both f/ < and V <Q 
on the U-V plane, in the simulation, are the star-forming 
gas in X < 5 kpc and F > 0. This area corresponds to 
the southern hemisphere for which observation is impos- 
sible with either the VLBA or VERA. 

4. DISCUSSION 

4.1. Rotation Speed of Star-Forming Regions 

iReid et al] ()2009bf ) investigated the peculiar motions of 
high mass star-forming regions in our Galaxy to obtain 
the Galactic parameters (i?o and Oq). They claimed that 
the star forming regions moves slower than the Galactic 
rotation. However, this is not the case in our result. As 
we can see in Figure [3 there is no clear tendency that 
the star-forming gas moves slower th an the true rotatin 
curve. ^° The observed slow rotation ()Reid et al.ll2009bl ) 
is probably caused by an uncertainty of the solar motion 

^'^ IReid et al.l l)2009bl ^ adopted the solar motion obtained by 
DB98 and that they introduced a model in which the star-forming 
regions move with a certain systematic velocity, K>fr, of —15 
km s~^with respect to the Galactic rotation. 
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Fig. 3. — : Top left: Distribution of stars projected onto the x-y plane (orange colors indicate stellar density). A major 
axis of the stellar bar inclines toward 25° from a y-axis (vertical axis). Top right: An edge-on view (x-z) of the stellar 
disk. Bottom left: Young stars, indicated by white-light blue color (age < 50 Myr) overlaid on old stars. Bottom 
right: Same as the left panel, but cold gas (white-green-black, Tg < 100 K) in stead of young stars. The unit length 
is in kpc. 



relative to the LSR, especially, Vq. If we compare with 
the right and left panels in Figure 1, it is clear that the so- 
la r motion of FA09 r educes the amount of slow rotation. 
If iReid et all (|2009b[ ) adopt the solar motion of FA09 in- 
stead of DB98, Vsfr could be smaller in their analysis. 
We have to note that "SUN" for an observer might have 
a large uncertainty in the peculiar motion. Therefore 
the peculiar motions of the star-forming regions could 
have a bi as with respect to th e Gala ctic rotation curve. 
Recently, iMcMillan fc BinnevI (|2009[ ) suggests that the 
value of Vq determined in these papers may be underes- 
timated by ~ 6 — 7 km s~^. Their value is close to the 



value of FA09 adopted in o ur paper. Our result is con - 
sistent with the argument o fMcMillan fc BinnevI (|2009l ). 

4.2. Origin of Large Peculiar Motions 

In order to discuss the origin of the large peculiar mo- 
tions of the star- forming gas, we compare three runs 
(Runs FID, ISO, and NSF): In Run FID, which is pre- 
sented in previous sections, the gas can cool down to 
20 K, and star formation from the dense and cold gas 
(riH > nth = 100 cm-3 and Tg < Tth = 100 K) is in- 
cluded. In Run ISO, the gas can cool down to 10"^ K, 
and star formation from the gas with njj > 0.1 cm"'^ 
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X [kpc] X [kpc] 

Fig. 4. — : Left: The peculiar velocities of randomly selected young stars (age < 10 Myr) represented by white arrows 
overlaid on distribution of stars projected onto the x — y plane (orange colors indicate stellar densities) at i = 1.25 Gyr. 
The peculiar velocity of a star particle is defined as the deviation from the true rotation curve. Right: The peculiar 
velocities of randomly selected star-forming gas (nn > 100 cm^^^ and Tg < 100 K) represented by red arrows are 
overlaid on distribution of cold gas (Tg < 100 K, green colors indicate their number densities). Note that very large 
peculiar velocities 100 km s^^) of a star-forming gas are caused by nearby supernovae explosions. Young stars and 
star- forming gas within i? = 3.0 kpc are not displayed. 



O 




5 10 15 

R [kpc] 

Fig. 5. — : Average rotation velocity of the star-forming 
ffas as a function of a galactocentric radius (solid line). 
Spacial distribution of these star-forming gas is shown in 
the right panel of Figure ID Each bin includes several 
tens samples, and error bars are the standard deviations 
of the mean. The dotted line indicates the true rotation 
curve. 



and T < 15000 K is included. In Run NSF, the gas can 
cool down to 20 K as in Run FID, but this run included 
neither star formation nor SNe. 

Figure [Hd shows the histogram of Wpcc/K:ir of the star- 
forming gas in Run ISO. They have peculiar velocities of 
~ 5% of the circular velocities, which cannot explain the 
observed large peculiar velocities of ~ 10% of the circular 
velocity. Figure [8^ shows the result of Run FID, where 
the star-forming gas has 2 times larger peculiar veloci- 
ties than those in Rim ISO. In order to understand what 
causes this difference, the distribution of a ratio of the 
pressure gradient (-Fkydro) to the non-axisymmetric com- 
ponent of the gravitational field (i^asym) was investigated 
as shown in Figure [9l In Run ISO, many star- forming 
gases are distributed around Fhydro ^ ^asym- However, 
most of the star-forming gases in Run FID are distributed 
-Fhydro < 0.1 X i^asym- Therefore, the non-circular motion 
of star-forming gas in Run ISO is not damped by hydro- 
dynamic pressure in the case of Run FID, and so, they 
can interact with transient stellar spirals via mechanism 
such as swing amplification (see i )3.1[) . 

Finally, we compare Run FID with Run NSF. The his- 
togram of Wpoc/K;ir of the star-formiug gas in Run NSF is 
shown in Figure [St. The distribution is almost the same 
as that in Run FID, suggesting that the acceleration by 
the SNe is not the main contributor in causing the large 
peculiar motions of the star-forming gas. 

4.3. Kinematic Distances and Face-on View of our 
Galaxy 
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Fig. 6. — : Left: Pseudo-observation of peculiar velocity of randomly selected star-forming gas (nn > 100 cm~^) 
relative to the "SUN" . The peculiar velocity is defined by the deviation from the flat rotation curve with a galactic 
rotation velocity at the "SUN". The position angle of the "SUN" is assumed to be 25° from the major axis of the 
stellar bar, and its galactocentric radius Rq = 8.0 kpc. The color represents density of cold (Tg < 100 K) gas. 
The star-forming gas within R = 3.0 kpc are not displayed. Right : The Peculiar velocities of the star-forming gas 
(riH > 100 cm~'^) in 5 kpc < i? < 10 kpc on a J7 — plane. The U and V are normalized to the galactic rotation 
velocity at the Sun in the model (triangles) and observations (open circles and crosses). For observational data, the 
solar motion relative to the LSR is the value of FA09 (open circles), and the value of DB98 (crosses). 



Since the pioneering work by I Port et all (|1958D there 
have been many attempts to construct the large-scale 
spiral structure of our Galaxy using the conventional 
"kinematic distance" method as a first approximation. 
The kinematic distance method is based on assump- 
tions that the orbits are not very different from circu- 
lar with an angular orbital speed decreasing monotoni- 
cally as a function of galactocentric radius. Nearby ma- 
jor spiral arms were identified by applying this method 
to v a,rious sources su c h as HI, CO, and HII r egions 
(e.g. lOort efaLl fl958|: iNakanishi fc Sofu3 [200l [20061 : 



iGeoreelin fc Georgelinl[l97l T 

The typical peculiar velocity of 20 — 30 km s~^, which 
is suggested by our numerical model and recent VLBI 
observations, might have a strong impact on how we see 
our Galaxy. The two panels of Figure [10] demonstrate 
the difference between the true gas distribution in our 
model galaxy and the view obtained using the kinematic 
distance technique, where the velocity field of gas is used 
to derive the distance from the "SUN" Since our aim 
here is to demonstrate how the large non-circular motions 
affect the estimations in kinematic distances, the errors 
caused by the distance ambiguity in the inner galaxy 
have to be in a minimu r n. Th erefore, we adopted the 
same method as iGomezj (|2006[ ) : based on the real dis- 
tribution of the ISM, we first figure out whether a tar- 
get gas element (particle) is located at the near or far 
side, and then we derive distance of the element using 
its line of sight velocity assuming the circular rotation. 



The reconstructed map roughly reproduce the real dis- 
tribution of the gas, but it shows many clear 'spur'-like 
structures towards the Sun, similar to those seen in the 

J ireviously published HI and CO niaps of our Galaxy 
Port et all 119581 : INakanishi fc Sofud [20031 . [200l . Most 
of these structures are not real, but they are the result of 
the large non-circular motions of the gas clouds as well 
as the non-axisymmetric flow caused by the bar. In ad- 
dition, outside the solar circle the arms in the left half of 
the figure have a pitch angle larger than that of the real 
arm, while in the right half the pitch angles in the kine- 
matic map are smaller than the real ones. This asym- 
metry is also visible in the maps constructed from the 
observations, suggesting that the motions of gas clouds 
in our Galaxy actually have large non-circular motions, 
similar to what was obtained in the present simulations, 
and that the location of the Sun relative to the major 
axis of central bar has an offset. More quantitatively. 
Figure [Tl] shows that the kinematic distances are well 
correlated with the real ones, but there are large errors 
with the typical value of ~ ±(2 — 3) kpc in the kinematic 
ones from the real ones. 

If we apply the kinematic distance method to the 
objects associated with the spiral arms, we must con- 
sider large non-circular motions. There is an approach 
which involves the modeling of non-circular motions 
of the gas instead of assum i ng a cir c ular motion (e.g. 



Foster fc Mac Will 
2008D. Recently, 



am; 



Gomed [20061 : iPohl et all 
Pohl et all (|2008D recovered a "real" 



1 [200 



ICombej 1 119911) showed a similar comparison bascd-on their 
'cloud-fluid' model, but a flxed stellar bar potential was assumed. 



3-D distribution of molecular gas from CP observations 
using a hydrodynamic simulation. They used a proba- 
bilistic method to match the observational CO gas dis- 
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Fig. 7. — : Same as FigurelHl but for (a) i?o = 7 kpc and (b) t = 1.35 Gyr. 



tribution to the model predictfon along the line-of-sight, 
instead of assuming a circular motion for the underly- 
ing kinematic model. H owever, the v e locity field of the 
ISM in previous works (|G6mezl l2006t iPohl et all l2008l) 
was probably much smoother than the real one because 
the ISM was treated as an isothermal gas (Tg ^ K) 
and the stellar bar was treated as a fixed potential. As 
discussed in H4:.2[ the cold gas clumps can have large non- 
circular motions through the non-linear gravitational in- 
teraction between them and transient stellar spirals. In 
order to recover the gas distribution in our Galaxy, it 
would be necessary to use a more complex velocity field 
as naturally expected for the multi-phase ISM in a live 
stellar disk. 



5. SUMMARY AND CONCLUSION 



High-resolution, A^-body-|-hydrodynamical simulation 
in which the multi-phase ISM and star formation are self- 
consistently taken into account revealed the kinematics 
of star-forming regions in a barred-spiral, like our Galaxy. 
We found that stellar and gas arms are not steady but 
arc transient and recurrent. This transient nature of the 
spiral arms may contribute to the large and complex non- 
circular motions of star-forming regions, as found in re- 
cent VLBI observations. 

We found that the peculiar motions of the Galactic star 
forming regions are qualitatively consistent with the mo- 
tions of the star forming gas in our model. Even though 
the kinematic distance technique roughly follows the real 
distance, each spiral arm is stretched as large as 4 — 6 
kpc along the line of sigt in the reproduced map. This is 
because of the assumption in the kinematic distance, i.e. 
the pure circular rotation. 
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Fig. 8. — : Distribution of the ratio of peculiar velocity 
to circular velocity, Wpcc/Kir, for star-forming gas in 5 
kpc < R < 10 kpc at t = 1.30 Gyr. The peculiar velocity 
(vpoc) is defined as the deviation from the local circular 
velocity, a) Run FID, b) Run ISO, and c) Run NSF. 
See text in M.'2\ for the detailed explanation about these 
Runs. 



The results shown here suggest that systematic as- 
trometric measurements of objects in the galactic disk 
are essential in order to understand the structure of our 
Galaxy. The VLBl astrometry and future space mis- 
sions for infrared astr ometry such as the original plan of 
Japanese JASMINE (jGouda et al.ll2006[ ) or even the as- 
trometry of nearby (a few kp cs) disk stars which sh ould 
be possible with the GAIA (jPerrvman et al.|[200l[) will 
give us a completely new picture of the spiral arms of our 
Galaxy, with the help of numerical simulations like the 
ones presented here, or ultimately cosmological galaxy 
formation simulations with sufficiently high enough res- 
olutions. 
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Fig. 9. — : Distribution of the ratio of hydrodynamical 
force (J^hydro) to asymmetric gravitational force (Fasym) 
for star-forming gas shown in Figure [H] Thick and thin 
curves indicate the results of Runs FID and ISO, respec- 
tively. 
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